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SUMMARY 


Constant amplitude tension-tension fatigue tests were 
conducted on AS4/3501-6 graphite epoxy (O2/02/-02)s laminates, 
where 0 was 15, 20, 25, or 30 degrees. Fatigue tests were conducted 
at a frequency of 5 Hz and an R-ratio of 0.1. Dye penetrant enhanced 
X-radiography was used to document the onset of matrix cracking in 
the central -0 degree plies, and the subsequent onset of local 
delaminations in the 0/-0 interface at the intersection of the matrix 
cracks and the free edge, as a function of the number of fatigue 
cycles. 

Two strain energy release rate solutions for local 
delamination from matrix cracks were derived: one for a local 
delamination growing from an angle ply matrix crack with a uniform 
delamination front across the laminate width, and one for a local 
delamination growing from an angle ply matrix crack with a 
triangular shaped delamination area that extended only partially into 
the laminate width from the free edge. Plots of G ma x vs. N were 
generated to assess the accuracy of these G solutions. The influence 
of residual thermal and moisture stresses on G were also quantified. 
However, a detailed analysis of the G components and a mixed-mode 
fatigue failure criterion for this material may be needed to predict 
the fatigue behavior of these laminates. 


INTRODUCTION 


In part I of this study [1], the damage observed in quasi-static 
tension tests on (O2/02/ _ 02)s AS4/3501-6 graphite/epoxy 

laminates, where 9 was 15, 20, 25, and 30 degrees, was documented. 
In addition, the stresses responsible for the matrix cracking and 
associated local delamination observed were calculated using finite 
element analyses. 

In part II of this study, fatigue tests were conducted on 

(O2/O2/-02)s AS4/3501-6 graphite/epoxy laminates, where 0 was 
15, 20, 25, and 30 degrees. Constant amplitude, tension-tension 
fatigue tests were conducted at a frequency of 5 Hz and an R-ratio 
of 0.1. Dye penetrant enhanced X-ray radiography was used to 

document the onset of matrix cracking, and the subsequent onset of 
local delaminations at the intersection of the matrix cracks and the 
free edge, as a function of the number of fatigue cycles. Two strain 
energy release rate, G, solutions were introduced for local 
delamination from matrix cracks. These G solutions were used to 

correlate the onset of local delamination in fatigue between layups 

with different off-axis angle plies. 

EXPERIMENTS 


Materials 

As documented in part I [1], panels were manufactured in an 
autoclave at NASA Langley Research Center from Hercules 
AS4/3501 -6 graphite epoxy prepreg using the standard 

manufacturer's procedure with a maximum temperature of 177°C 
(350°F) in the cure cycle. One panel was produced for each of the 
(O 2 /02/-02)s layups, where 0 was 15, 20, 25, or 30 degrees. Panels 
305 mm (12 in.) square were manufactured and then cut into 20 
individual coupons, 127 mm (5.0 in.) long by 25 mm (1.0 in.) wide. 
Each specimen had a nominal ply thickness of 0.124 mm (0.0049 in.) 
and a fiber volume fraction of 67%. 

Experimental Procedures 

All static tests were conducted using the procedure outlined in 
part I of this study [1]. Fatigue specimens were placed in the 
hydraulic grips of a hydraulic load frame such that each specimen 
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had a gage length between the grips of approximately 76 mm (3.0 
inches). A 25 mm (one inch) gage length extensometer was mounted 
on the specimen midway between the grips to measure nominal 
strain. As specimens were loaded quasi-statically to the mean load, 
a plot of load versus nominal strain was generated using an X-Y 
recorder. The extensometer was then removed, and the specimens 
were loaded at a cyclic frequency of 5Hz and an R-ratio of 0.1. Tests 
were run at several maximum load levels below the critical load 
required to initiate damage in the static tests. The cyclic loading 
was stopped at fixed cyclic intervals and a zinc iodide solution that 
was opaque to X-rays was placed on the edge using a hypodermic 
needle. A radiograph was taken of the specimen width, with the 
specimen loaded at the mean load, to document the presence of 
matrix cracking and delamination. 

Table 1 lists the maximum cyclic loads and the schedule used 
to obtain radiographs for each cyclic load level. For example, tests 
run at maximum cyclic loads of 5000 and 6000 pounds were 

interrupted every 1000 cycles to obtain a radiograph of the 
specimen. Tests run at a maximum cyclic load of 3500 lbs were 
interrupted every 1000 cycles for the first 10,000 cycles; every 

2,000 cycles between 10,000 and 20,000 cycles; every 5,000 cycles 
between 20,000 and 50,000 cycles; and every 10,000 cycles 

thereafter, to obtain a radiograph of the specimen. Loading was 
continued until a radiograph clearly indicated the presence of a local 
delamination associated with a matrix crack near the free edge. 
After the fatigue specimens were removed from the hydraulic grips, 
the edges were polished and examined in the light microscope. 

Photographs of the polished edges were taken through an optical 
microscope to document the damage observed on the straight edge. 

Initially, some coupons were coated with a thin film of water- 
based typewriter correction fluid on either edge. After a few fatigue 
tests were performed, however, it was apparent that the onset of 
matrix cracking and local delamination was difficult to detect 
visually on the edge with this brittle coating. Hence, most of the 
fatigue coupons were tested without coating the edges. It was also 
easier to see the formation of the matrix cracks in the radiographs 
if there was no coating on the edge, because the dye penetrant used 
to reveal the damage was more likely to accumulate near a coated 
edge. 
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Experimental Results 
Static Tests 


The results from the quasi-static tests were documented in 
part I of this study [1], The radiographs showed triangular-shaped 
shaded regions bounded by the free edge and matrix cracks, 
indicating that a local delamination was present. Edge micrographs 
showed that local delaminations were present in both 0/-9 
interfaces, extending in the same direction from curved matrix 
cracks in the central -0 degree plies. 

Fatigue Tesfs 

The same type of damage seen in the radiographs and edge 
micrographs of specimens taken after the quasi-static tests was 
present in the specimens that were tested in fatigue [1], However, 
unlike the quasi-static tests, the sequence of damage was evident in 
the series of radiographs taken during the cyclic loading. For 
example, fig. 1 shows radiographs taken immediately after the onset 
of matrix cracking (fig. la), and the subsequent onset of local 
delamination from the matrix crack (fig. 1b), for a (0 2 /20 2 /-20 2 ) s 
laminate that was subjected to a maximum cyclic load of 15,570 N 
(3500 lbs). Similar radiographs were obtained for the other layups 
tested in fatigue. At the number of cycles corresponding to local 
delamination onset, all of the specimens had at least one, and in 
most cases, several matrix cracks in the central -0 degree plies. A 
few matrix cracks were also observed in the +0 degree plies for 
some laminates. Table 2 lists the average number of matrix cracks 
in the gage section measured from the radiographs of laminates 
taken just after the onset of local delamination. Local delaminations 
were present in the 0/-0 interfaces near some, but not all, of these 
matrix cracks. By the time the fatigue tests were terminated, 
matrix cracks had also developed in the +0 degree plies of many 
laminates. Delaminations in the 0/-0 interfaces often terminated at 
these +0 degree matrix ply cracks (fig. 2), 

Figures 3-6 show the cumulative number of cycles to the onset 
of matrix cracking and subsequent local delamination onset for each 
layup. The load at onset of damage (matrix cracking and local 
delamination) under a monotonically increasing load, determined in 
part I of this study, is shown on the ordinate. Figures 7-10 show the 
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cumulative number of cycles to local delamination onset minus the 
number of cycles to onset of matrix cracking. Hence, figures 7-10 
represent the number of cycles to delamination onset once the 
matrix crack is present. 

STRAIN ENERGY RELEASE RATE ANALYSIS 


Uniform Local Delamination 


In reference [2], an equation was derived for the strain energy 
release rate associated with a uniform local delamination growing 
from a 90 degree matrix crack (fig. 11a) in a multidirectional 
laminate using the classical definition for G as 


2 dA 


(1 ) 


where dC is the incremental change in laminate compliance 
associated with an incremental change in surface area, dA, created 
by the local delamination. This solution assumed that the composite 
consisted of a laminated region and a locally delaminated region in 
series, resulting in the following expression 


2mw 2 


1 L_ 

t|_DELD tE LA M 


where P is the applied mechanical load on the laminate, w is the 
laminate width, and m is the number of delaminations growing from 
the matrix crack. The parameter, m, would have a value of 2 if the 
cracked off-axis ply is in the interior of the laminate, corresponding 
to local delamination on either side of the matrix crack. However, m 
would have a value of 1 if the cracked off-axis ply was a surface 
ply. In eq(2), t is the laminate thickness, Elam is the laminate 
modulus (as calculated from laminated plate theory), tLD is the 
thickness of the locally delaminated region that carries the load 
(i.e., the laminate thickness minus the thickness of the cracked off- 
axis plies), and Eld is the modulus of the locally delaminated region 
(as calculated from laminated plate theory). 

As shown in appendix 1, eq(2) also applies to the case of a 
uniform delamination growing from an angle ply matrix crack as 
shown in fig. (11b). For the laminates in this study, t = Nh and tLD = 
(N-n)h, where N is the number of plies in the laminate, n is the 
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number of cracked off-axis plies, and h is the ply thickness. 
Therefore, noting that N x = P/w, eq(2) may be rewritten as 


2mh \(N-n)E|_D NElam; ’ 1 

The twelve-ply (O2/02/-02)s laminates tested in this study may be 
analyzed as six-ply laminates with twice the measured ply 
thickness. Hence, N=6, n=2, and N-n=4. Then, noting that m=2 for an 
interior -0 matrix ply crack, a normalized strain energy release rate 
may be calculated as 


a.ilJ i_| 

isj2 4 UEld 6E L am/ 


(4) 


Bending-Extension Coupling Influence on G 

The value of Eld in eq.(4) depends on whether or not the locally 
delaminated region consists of symmetric or asymmetric 
sublaminates. For example, the locally delaminated region may be 
thought of as a (0/0 ) s symmetric sublaminate. In this case, the 
modulus of the locally delaminated region is simply 


Eld = 


1 

(N-nMAn )' 1 


(5) 


where (An)- 1 is the first term of the inverse extensional stiffness 
matrix for the (O/0) s sublaminate. However, as noted in ref. [3], if 
the locally delaminated region consists of two 2-ply asymmetric 
sublaminates, then the presence of bending-extension coupling 
should be reflected in the modulus calculation. This may be 
accomplished using laminated plate theory by prescribing an applied 
mechanical strain, e x , setting the applied N y , N xy , and M y equal to 0, 
and setting k x and K xy equal to 0. The N x that results from this 
analysis is derived in appendix 2 and may be used to determine Eld 
as 


6 



( 6 ) 


Eld = 


2 

(N-n)h 


An + 



where Di, D 2 , D 3 , and D are functions of various terms in the A, B, 
and D matrices as defined in appendix 2. 

Any difference between equations 5 and 6 reflects the 
constraint imposed on the 2-ply asymmetric sublaminate by the 
undamaged laminate and the grips (if the sublaminate extends into 
the grips). Because of this constraint, the sublaminate is not able to 
develop a k x or K xy curvature. However, under tension loading, a non- 
zero K y curvature may develop away from the grips in the 
sublaminate. For the type of sublaminates that form as a result of 
edge delamination [3-5], this influence of bending-extension 
coupling may be significant. However, for the (0/6 )t sublaminates 
that constitute the locally delaminated region in the (O/9/-0 ) s 
laminates with -0 matrix cracks, the influence of bending-extension 
coupling is small. This small influence is shown in fig. 12, where the 
normalized G of eq(4) is calculated as a function of 0 for (O/0/-0) s 
laminates with -9 matrix cracks using material properties in Table 
3. The normalized G values determined by calculating Eld from eq(5) 
assuming a symmetric (O/0) s sublaminate are nearly identical to the 
normalized G values determined by calculating Eld from eq(6) using 
the loading prescribed above in the laminated plate theory analysis. 
However, unlike many of the sublaminates formed by edge 
delamination, the (O/0)j sublaminates that constitute the locally 
delaminated region in the (O/0/-0) s laminates with -0 matrix cracks 
will also exhibit shear-extension coupling. This shear-extension 
coupling may have more influence on the calculated Eld than the 
bending extension coupling. 

Shear- Extens/on Coupling Influence on Q 


The influence of shear-extension coupling on the calculated 
Eld for the 2-ply asymmetric sublaminate may be determined using 
laminated plate theory by prescribing the loading outlined above, but 
setting y xy = 0 instead of prescribing N xy = 0. The N x that results 
from this analysis is derived in appendix 2 and may be used to 
determine Eld as 


Eld = 


2 

(N-n)h 



/ A12D22 ~ B12B22 
\A22D22 - B22B22 


A 12 + 


( A 12 B 22 - A 22 BM B 
\A 22 D 22 - B 22 B 22 ' 


(7) 
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This Eld value reflects the shear constraint imposed on the locally 
delaminated region by the undamaged laminated region. The shear 
constraint results in a greater apparent stiffness, Eld, and hence, a 
correspondingly lower strain energy release rate. Figure 12 shows 
the normalized G of eq(4), calculated as a function of 0 for (O/0/-0) s 
laminates with -0 matrix cracks, where the normalized G values 
were determined by calculating Eld using eq(7). This loading 
assumes that the locally delaminated region consists of two 2-ply 
(O/0 )t asymmetric sublaminates with bending-extension and shear- 
extension coupling. As shown in fig. 12, the normalized G values 
determined by calculating Eld using eq(7) are significantly lower 
than the normalized G values determined by calculating Eld using 
eqs(5&6) for 0 between 15 and 30 degrees. 

Uniform Locdl Dela m instion vs. Edge Delamination 

As noted earlier, the delaminations that form in the 0/-0 
interfaces of (O/0/-0) s laminates are localized and bounded by -0 
degree ply matrix cracks. Hence, the physical evidence dictates that 
these delaminations should be analyzed as local delaminations. 
However, it may be enlightening to calculate the strain energy 
release rate, G, for edge delamination in the 0/-0 interfaces of the 
(O/0/-0) s laminates and compare that to the corresponding G for 
local delamination. The strain energy release rate, G, for edge 
delamination is given by [4] 


G (Elam - E ) (8) 

where e is the nominal applied strain, t is the laminate thickness, 
Elam is the laminate modulus, and E* is the modulus of the laminate 
once it is totally delaminated, as determined by the rule of mixtures 


M 

I Ejtj 

E*=-^!-_ (9 , 

where M is the number of sublaminates formed by the edge 
delamination and Ej and tj are the modulus and thickness of the i 
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sublaminates, respectively. For edge delamination in the 0/-0 
interfaces of the (O/0/-0) s laminates, eq(9) becomes 

E ._2E Hi jH S 

3 

The strain energy release rate calculated in eq(8) will depend on the 
assumptions used to determine the (0/0) and (-0) sublaminate moduli 
in eq(10), which determines the delaminated modulus, E‘, in eq(9). 

To compare the different possibilities, eq(8) was rearranged 
to yield a normalized G as 


— = -^Elam - E ) 
e 2 t 2 


Figure 13 shows the normalized strain energy release rate, G/e 2 t, 
for edge delamination in the 0/-0 interfaces of the (O/0/-0) s 
laminates as a function of 0 where the sublaminate moduli were 
determined assuming they were either symmetric, asymmetric with 
bending-extension coupling only (B/E Coupling), or asymmetric with 
shear-extension coupling in addition to bending-extension coupling 
(B/E & S/E Coupling). Unlike the comparison for local delamination, 
there is a significant difference in the normalized G for edge 

delamination between the symmetric sublaminate case and the 

bending-extension coupling case. This difference results from the 
effect of bending-extension coupling on the modulus of the -0 
sublaminate in eq(10) which does not contribute to G for local 

delamination. Furthermore, when the influence of shear-extension 
coupling is included, the normalized G for edge delamination 

becomes very small for all values of 0. 

In order to compare normalized G values for both local and 
edge delamination, Hooke’s law in the form N x = NheE|_AM , where N=6, 
was substituted into eq(4) to yield 


3E L am j 1 

_ 1 ] 

2 Ue ld 

6E|_ami 


Figure 14 shows G normalized by e 2 t for edge delamination and 
uniform local delamination, as calculated by eqs.(11) and (12) 
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respectively, as a function of 0 for delamination in the 0/-0 
interfaces of the (O/0/-0) s laminates. The sublaminate moduli were 
determined assuming they were asymmetric, with shear-extension 
coupling in addition to bending-extension coupling. Figure 14 
indicates that, for the same applied strain, the strain energy release 
rate associated with a uniform local delamination from a matrix 
crack exceeds the strain energy release rate associated with edge 
delamination for all values of 0 . 

Partial Local Delamination 

The local delamination growing from an angle ply matrix crack 
modeled in fig ,11b has a uniform delamination front across the 
laminate width. However, the local delaminations observed in the 
experiments have triangular shaped delamination areas that extend 
only partially into the laminate width from the free edge [1], 
Therefore, an analysis was developed to account for the fact that 
the local delaminations in the (O/0/-0 ) s laminates did not extend 
across the entire laminate width, but instead were limited to the 
extent of the matrix crack in the -0 degree ply. The model was based 
on the assumption that the regions near the free edges containing 
the partial local delaminations act as a system of sublaminates in 
series, just as was assumed for the uniform local delamination, but 
this series system was assumed to be in parallel with the 
undamaged laminate in the interior of the laminate width (fig. 11c). 

Figure 15 shows three possible ways of visualizing the 
delamination growing from the matrix crack in the -©degree ply. Fig. 
15a shows a triangular shaped delamination that intersects the free 
edge at 90 degrees. Fig. 1 5b shows a triangular shaped delamination 
that intersects the matrix crack at 90 degrees. Finally, fig. 15c 
shows a triangular shaped delamination that intersects the free 
edge at an angle of 20. Eventually, if the local delamination in the 
0/-0 interface reaches a +0 degree matrix crack and arrests, the 
local delamination will appear to be an isosceles triangle with +0 
and -0 angles intersecting the straight free edge. The relationship 
between the matrix crack length, d, and the delamination length 
along the free edge, a, is different for the three cases shown in fig. 
15. These relationships are given in Table 4. Hence, the shape of the 
delamination front must be known in order to relate the 

delamination length on the free edge to the length of the matrix 
crack. 



For a local delamination from a matrix crack, the compliance, 
S, where S = 1/E, for a laminate with a given local delamination 
area, A, in an interface whose total area is A*, may be represented 
as (appendix 1) 

S - Slam _ A 

S* - Slam A* (13) 


where S‘ - 1 /E*ld . and E'ld = (tLDEi_D)/t as defined in ref. 2. 
Assuming a partial local delamination is present at a matrix crack 
on both free edges (fig. 11c), and noting that a'=ka, where k is given 
in Table 4 for the three geometries shown in figures 15 a-c, then 
A = maa’ = mka 2 and A* = 2mLka. Substituting for A and A* in eq(13) 
for any of these three geometries yields 


s = 1/ 1 

2uldEld 



1 

Elam 


(14) 


Recalling from ref. 4 that, for a system in parallel, 

G = - J- Ve 2 
2 dA 


(15) 


and noting that for the partial local delaminations shown in fig. 11c 
A = mka 2 , where k is given in Table 4 for the three geometries 
shown in figures 15 a-c, then dA = 2mkada. Substituting V = wLt and 
dA into eq(15) yields 


G = -Iwid(L)dE 

2 2mk da 


(16) 


The modulus of a laminate containing a delaminated region near the 
free edges loaded in parallel with an undamaged region in the 
interior of the laminate width is given by [4] as 

E = (E* - Elam) + Elam 

(a) (17) 

For the partial delamination case shown in fig. 11c, A' = 2ml_a' and 
(A’)* = mwL. Substituting these expressions in eq(17), and noting 
that a'=ka, where k is given in Table 4 for the three geometries in 
fig. 15a-c, yields 


(18) 


E - (e* - Elam) ^ a + Elam 


Unlike E* calculated for edge delamination in eq(9), E* in eq(18) is a 
function of the local delamination length along the free edge, a. 
Differentiating eq(18) with respect to a, and substituting into 
eq(16) yields 



(19) 


Substituting E = 1/S, where S is given by eq(14), into eq(19) and 
differentiating yields 


where 


and 



— 1 + Ba 

(Ba + C) (Ba + C) 2 


( 20 ) 


B =lf — 1 — 
2 UldEld 


1 

iElam 



(21) 


C = 


1 

Elam 


( 22 ) 


Although the relationship between the matrix crack length, d, and 
the delamination length along the free edge, a, is different for the 
three cases shown in fig. 15, eq.(20) is independent of k, and hence, 
G(a/I_) is independent of the local delamination geometry assumed. 
However, as noted previously for eq(2), G in eq(20) will depend on 
the assumptions made when determining Eld in eq(21). Therefore, B 
in eq(21), and hence G in eq(20), will depend on whether eqs. 5,6, or 
7 are used to calculate Eld- For the remainder of this paper, it will 
be assumed that the sublaminates that constitute the locally 
delaminated region experience both bending-extension and shear- 
extension coupling. Hence, eq(7) will be used to calculate Eld. 

Eq.(20) may be rearranged to yield 


G = p- 

2m 



1 2 


(23) 


where D = BL. When a/L = 0, then eq.(20) becomes 


G = 


^L/2D\ 

2m \ C 2 j 


( 24 ) 


Substituting for D and C, normalizing by e 2 t, and noting that ti_D = 4h 
and t=6h, eq(24) becomes 


_G_ = 3E lam /_! 1 | 

e 2 t 2 UEld 6Elam/ (25) 

which is identical to eq(12) for the uniform local delamination. 

Figs. 16-19 compare normalized G values calculated using 
eq(23) for partial local delamination, as a function of delamination 
length along the free edge normalized by the laminate length, a/L, 
to normalized G values calculated from eq(12) for uniform local 
delamination and from eq(11) for edge delamination modeled in the 
0/-0 interfaces of (O/0/-0) s laminates, where 0 = 15, 20, 25, and 30 
degrees, respectively. The uniform local delamination and edge 
delamination solutions are independent of delamination length, but 
for the partial local delamination solution, G decreases with 
increasing delamination length along the free edge. Both local 
delamination solutions yield greater normalized G values than the 
edge delamination solution for all delamination lengths. 

Fig .20 shows G normalized by e 2 t, as a function of a/L up to 
a/L=0.1, calculated using eq(20) for partial local delaminations in 
the 0/-0 interfaces of (O/0/-0) s laminates, where 0 = 15, 20, 25, and 
30 degrees. Normalized G values were only plotted for a/L values up 
to 0.1 because the magnitudes for small matrix crack lengths, and 
hence small delamination lengths, were of primary interest for 
comparing to experimental data. In each case, G decreases slightly 
with delamination length, unlike the solutions for uniform local 
delamination which were independent of delamination length. If the 
partial delamination geometry can be anticipated and linked to the 
matrix crack length as shown in fig.15a-c, then both the laminate 
toughness as a function of fatigue cycles and a critical matrix crack 
length would be needed to predict the onset of a partial local 
delamination in fatigue. However, as shown in fig. 20, if the critical 
matrix crack length is sufficiently small, and hence, the 
corresponding delamination length along the free edge is small, the 
difference between the uniform and partial delamination solutions 
will be insignificant. In this case, G may be calculated from eq(2) to 
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compare to experimental data, even though the local delamination is 
not uniform across the laminate width. 

Thermal/Moisture Contribution to G 

The strain energy release rate for local delamination will be 
influenced by the presence of residual thermal and moisture 
stresses in the laminate. A uniform local delamination G solution 
that includes the effect of the mechanical load, the residual thermal 
stresses that exist in the laminate after it cools down from the 
curing temperature, and the hygroscopic stresses that exist in the 
laminate due to the presence of absorbed moisture, was recently 
developed [6]. Figs. 21 -24 show the contribution of mechanical load 
(M), temperature (T), and hygroscopic (H) moisture weight gain to G 
for a uniform local delamination in the e/-0 interfaces of (O/0/-6) s 
laminates (where 6 = 15, 20, 25, and 30 degrees, respectively) 
subjected to a mechanical load of 3000 Ib/in. and a AT of -156°C 
(-280°F). As shown in figs. 21 -24, when the contribution of the 
residual thermal stress is included, G^+T will increase above the G^ 
value obtained from equation (2) for an applied mechanical load 
alone. However, when the presence of absorbed moisture is included, 
as may occur from exposure to the ambient air (appendix 3), the 
thermal stresses will relax and reduce the strain energy release 
rate, G M+T+H , with increasing AH. Similar behavior was observed for 
edge delamination [7]. However, the magnitude of the thermal and 
moisture effects is much more significant for edge delamination. 
The laminates tested in this study had a moisture weight gain, AH, 
between 0.2 and 0.3 % (appendix 3). For the (0/e/-6) s laminates 
shown in figs. 21-24, this corresponds to a 5.8-7.7% increase in 
qM+t+h compared to the G M solution for mechanical loading only. 
The difference in G M+T+H and G M is even less at higher mechanical 
load levels (Table 5). 

CORRELATION OF ANALYSIS AND EXPERIMENTS 

Ideally, the closed form solutions for local delaminations from 
matrix cracks should be verified and calibrated by comparing them 
to detailed finite element solutions that sum the G components for 
the three fracture modes to obtain a total G. Such comparisons have 
been made for edge delamination [3-5]. However, to date only the 
case of uniform local delamination from 90 degree matrix cracks 
has been verified [8]. Part of the difficulty is that the modelling for 
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local delamination requires a full 3D analysis, with the matrix crack 
geometry modeled accurately. The solutions for non-90 degree 
matrix cracks that have been attempted to date [9,10], have assumed 
that the angle ply matrix cracks may be represented as straight 
cracks through the laminate thickness that are normal to the 02/-02 
interface. However, as seen in micrographs of the edges for the 
(O2/02/-02 )s laminates [1], these angle ply matrix cracks are curved, 
intersecting the ply interface at oblique angles. This difference in 
actual versus assumed matrix crack geometry through the laminate 
thickness may have a significant effect on the G components 
calculated from a finite element analysis for local delamination. 

Although the fatigue behavior of the (O2/02/-02)s laminates 
cannot currently be predicted, it may be possible to correlate the 
fatigue delamination onset behavior between laminates with 
different values of 0. The data from each layup may be used to 
generate delamination onset criteria by generating plots of 
maximum cyclic G versus the number of cycles to delamination onset 
for the composite material. A similar characterization has been used 
previously for edge delamination data [5,11]. 

Laminate and sublaminate moduli for the (O 2 / 0 2 /-© 2 ) s 
AS4/3501-6 graphite epoxy laminates were calculated using the 
properties in table 3 in laminated plate theory. The moduli of the 
locally delaminated regions, Eld. were calculated using eq(7) to 
incorporate the influence of bending-extension and shear-extension 
coupling. These moduli were then used in eq(2), along with the 
maximum cyclic loads in the fatigue tests, to calculate strain 
energy release rates for uniform local delaminations corresponding 
to the measured number of cycles to delamination onset for each 
layup. In addition, the matrix crack length, d, corresponding to the 
number of cycles for local delamination onset was measured from 
the X-ray radiographs where local delamination was first observed 
for each specimen. Examination of the partial local delaminations in 
these radiographs indicated that the geometry shown in fig. 15c was 
reasonable for the layups tested. Therefore, measured matrix crack 
lengths, d, were used in the expression in Table 4 corresponding to 
fig. 15c to determine critical delamination lengths along the free 
edge, a, at local delamination onset for each layup. Table 6 shows 
the range of delamination lengths along the free edge, a, normalized 
by the gage length, L, measured just after local delamination onset 
in fatigue. These critical delamination lengths were used, along with 
the maximum cyclic loads in the fatigue tests, to calculate strain 
energy release rates from eq(20) for partial local delaminations 
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corresponding for the measured number of cycles to delamination 
onset for each layup. Hence, plots of G m ax versus the number of 
cycles to delamination onset, N, were generated for each layup 
assuming either uniform, or partial, local delamination onset from a 
-0 degree matrix ply crack. 

Figure 25 shows the G ma x vs. N plot for both the uniform and 
partial local delamination models. The maximum cyclic G levels for 
the uniform and partial local delamination models are very similar 
because the experimentally determined critical delamination lengths 
were small (Table 6, fig. 20). In addition, for any particular 
maximum cyclic load, the maximum cyclic G levels for all values of 
0 also agree fairly well, indicating that the assumptions of shear- 
extension and bending-extension coupling reflected in eq(7) was 
accurate (fig. 12). Finally, the number of cycles to delamination 
onset at each cyclic G level appear to be greater for the 15 and 20 
degree cases than for the 25 and 30 degree cases. These differences 
in the number of cycles to delamination onset may be due to a 
difference in either the mixed-mode ratio, or the extent of matrix 
cracking prior to delamination onset, for the individual cases. 
Neither of these effects are incorporated in the current analyses. 

Figure 26 shows the G ma x vs. N plot for the uniform local 
delamination models with mechanical loading only, G M , and with 
mechanical, thermal, and hygroscopic loading combined, G M+T+H . For 
any particular maximum cyclic load, the maximum cyclic G m+t +h is 
only slightly greater than the the maximum cyclic G M because of the 
relatively small contribution of the residual thermal and moisture 
stresses to G (Table 5). 

The G max vs. N plots in figures 25 and 26 provide some insight 
into the accuracy of the G solutions derived in this paper. However, 
as noted earlier, a detailed analysis of the G components and a 
mixed-mode fatigue failure criterion for this material may be 
needed to predict the fatigue behavior of these (O 2 / 0 ?/- 0 p)c 
laminates. 


CONCLUSIONS 


Constant amplitude tension-tension fatigue tests were 
conducted on AS4/3501-6 graphite epoxy (O2/02/-02)s laminates, 
where 0 was 15, 20, 25, and 30 degrees. Fatigue tests were 
conducted at a frequency of 5 Hz and an R-ratio of 0.1. Dye penetrant 
enhanced X-ray radiography was used to document the onset of 
matrix cracking in the central -0 degree plies, and the subsequent 
onset of local delaminations in the 0/-0 interface at the 
intersection of the matrix cracks and the free edge, as a function of 
the number of fatigue cycles. Unlike the quasi-static tests [1], the 
sequence of damage was evident in the series of radiographs taken 
during the cyclic loading, with the matrix cracking always preceding 
the onset of local delamination. 

Two strain energy release rate solutions for local 
delamination from matrix cracks were derived: one for a local 
delamination growing from an angle ply matrix crack with a uniform 
delamination front across the laminate width, and one for a local 
delamination growing from an angle ply matrix crack with a 
triangular shaped delamination area that extended only partially into 
the laminate width from the free edge. Plots of G ma x vs. N were 
generated to assess the accuracy of these G solutions. The influence 
of residual thermal and moisture stresses on G were also quantified. 
However, a detailed analysis of the G components and a mixed-mode 
fatigue failure criterion for this material may be needed to predict 
the fatigue behavior of these laminates. 
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APPENDIX 1 - UNIFORM LOCAL DELAMINATION FROM AN ANGLE PLY 

MATRIX CRACK 

For a local delamination from a matrix crack, the laminate 
compliance, S, where S = 1/E, for a laminate with a given local 
delamination area, A, in an interface whose total area is A*, may be 
represented as 


S - Slam _ A 
S* - Slam A* 


( A 1 ) 


where S - 1/E ld . and E ld = (ti_DELD)/t as defined in ref.2. Fig. 11a 
shows the uniform local delamination that extends across the entire 
laminate width and grows normal to a 90 degree matrix ply crack in 
the load direction. For the uniform delamination in fig. 11a, A = 
mwa, and A* = mwL, where m=2 for an interior matrix crack and m=1 
for a matrix crack in a surface ply. Substituting these expressions 
into eq(A1) yields 


S = 


^ldEld 



1 

Elam 


(A2) 


which is identical to eq(7) in ref.2 for local delamination from a 90 
degree matrix ply crack. 

Fig. 11b shows the uniform local delamination that extends 
across the entire laminate width and grows from a -9 degree matrix 
ply crack normal to the load direction. For the uniform delamination 
in fig. 11b, A= (mw/sin9) (a sin0) = mwa, and A* = mwL. Substituting 
these expressions into eq(A1) also yields eq(A2). 

For the 90 degree matrix crack induced local delamination in 
fig 11a, dA = mwda. Similarly, for the -0 degree matrix crack 
induced local delamination in fig. 11b, dA = (mw/sin0) d(a sin0) = 
mwda. In ref.2, G is defined as 


G = 1 v ° 2 tt 

2 dA 


(A3) 


Substituting for V = wLt and dA = mwda and differentiating eq (A2) 
with respect to "a" yields 
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G =— E2— ( — I 1 — I 

2rnw 2 UldElD tLAMEL\M/ 


(A4) 


for uniform local delaminations from either a 90 degree or a -0 
degree matrix crack. 


APPENDIX 2 - ANALYSIS OF MODULI FOR ASYMMETRIC SUBLAMINATES 


In laminated plate theory, the load and moment resultants are 
related to the midplane strains and curvatures as 



(A5) 


where A, D, and B represent the extensional, bending, and bending- 
extension coupling matrices, respectively. By prescribing an applied 
mechanical strain, e x = eo, setting the applied N y , N xy , and M y equal to 
0, and setting k x and K xy equal to 0, eq(A5) becomes 


/ Nx ) 



£o 1 

0 



£y 

0 


A B 

1 

Yxy 

M x 

f 

B D 

1 

0 

0 



K y 

M x y j 



o I 


(A6) 


The three equations whose left hand side are zero may be written as 


-A]2£o = A22£y + A26Yxy + 1*22 K y 

■A]6E{) = A26^y + Afi^Yxy + I?26Ky 
-B ]2 e 0 = B 2 2£y + B 2 6Yxy + D 2 2 Ky 


(A7) 


where Ajj=Ajj. This set of three non-homogeneous equations may be 
solved for e y , Yxy. and K xy using Cramer's rule. Hence, 


Di D? 


= D3 
y D 


(A8) 
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where 


D,= 


■A 12 
"A j6 

■B 1 2 


A 26 
A 66 

B 26 


B 22 

B 26 

D 22 


D 2 = 


A22 'A 12 B22 
A 2 6 'A 16 B26 
B22 'B12 D22 



A 22 

A 2 6 

B22 


A26 -A 12 
A 66 "A 16 
B 2 6 -B )2 


and 



A 22 A26 B22 
A26 A 66 B26 
B22 B26 D22 


These expressions for ey, y xy , and K X y may be substituted into the 
first equation of eqs(A6) to yield 


N x = 


a,i+ (^) a '2 + (^) a '‘ + (§) b 


12 


eo 


(A9) 


Hence, N x may be obtained explicitly as a function of the terms in 

the A, B, and D matrices and the prescribed axial strain, eg Then 
noting that 


Eld = 


2N X 

(N-n)he 0 


( A 1 0) 


yields 


Eld = 


2 

(N-n)h 




(All) 


By prescribing the loading outlined above, but setting y xy = 0 
instead of prescribing thatN X y= 0, i.e., by prescribing an applied 
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mechanical strain, e x = eo, setting the applied N y and M y equal to 0, 
and setting y xy , k x , and K xy equal to 0, eq(A5) becomes 


/ N x | 



1 eo 

0 



£y 

N X y' 


A Bl 

0 

M x 

f “ 

.B DT 

0 

0 



K y 

> M xy 1 



0 / 


( A 1 2) 


The two equations whose left hand side are zero are 


0 = A]2£o + A22£y + B22Ky 


( A 1 3) 


0 = Bi2£o + B 22 £y + D22% 


( A 1 4) 


where Ajj=Ajj and By=Bjj. Multiplying eq(A13) by -D 22 . multiplying 
eq(A14) by B 22 . and then adding the two equations yields 

e -( B 12 B 22 - A 12 D 22 U 

y \A 22 D 22 - B 22 B 22 / (A15) 

Furthermore, multiplying eq(A13) by -B 22 , multiplying eq(A14) by 
A 22 . and adding the two equations yields 

K _ [ A 12 B 22 - A 22 B 12 I 

y U 22 D 22 - B 22 B 22 / (A16) 


Finally, substituting eqs(A15) and (A16) into the first equation of 
eq(A12) yields 


N x = 



Al2D22 ~ Bl2B22 | A , ( Al2B22 - A22Bl2 | g 
A22D22 ‘ B22B22' 'A 22 D 22 - B22B22/ 


eo 


( A 1 7) 


Hence, N x may be obtained explicitly as a function of the terms in 
the A, B, and D matrices and the prescribed axial strain, eo. Then 
recalling eq(A10) 


Eld = 


2N X 

(N-n)he 0 


and substituting the expression for N x from eq(A17) yields 
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Eld = 


(N-n)h 


^11 


/ A 12 D 2 2 - B12B22 ] A , ( A12B22 - A 2 2Bi2 \ g 
\A 22 D 22 - B 22 B 22 * \ A 2 2 D 22 - B 22 B 22 / 


(A18) 


APPENDIX 3 - MOISTURE CONTENT DETERMINATION 

Just prior to testing the (O2/02/-02)s laminates, two control 
coupons of each layup were weighed, dried, and then weighed again 
to determine AH for the laminates in the ambient environment. These 
dried control specimens were then weighed periodically to identify 
their moisture content by weight as a function of time. 

Fig. 27 shows the percentage moisture content by weight, AH, 
as a function of the number of days from the time the specimens 
were dried. There was very little difference in the percentage 
moisture weight gain for the different layups. The percentage 
difference in the original weight measured just before drying and 
the weight measured immediately after the specimens were dried 
was approximately 0.2% for all the layups. This value is shown in 
fig. 27 as the original ambient condition. Because all of the fatigue 
data were generated within a period of 30 days, the percentage 
moisture weight gain for all the laminates that were tested in 
fatigue was assumed to be between 0.2% and 0.3%. 
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TABLE 1 


X-ray Schedule for (O2/02/-02)s Fatigue Tests 


Pmax. N (lbs) 

X-ray taken 

Between Ni 



every N cycles 

cycles 


26,690 (6000) 

1000 

1 

5000 

22,240 (5000) 

1000 

1 

5000 

17,790 (4000) 

1000 

1 

10,000 


2000 

10,000 

20,000 

15,570 (3500) 

1000 

1 

10,000 


2000 

10,000 

20,000 


5000 

20,000 

50,000 


10,000 

50,000 

1 10,000 

13,340 (3000) 

10,000 

1 

100,000 


50,000 

100,000 

1,000,000 
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TABLE 2 - Average Number of Martrix Cracks at Local 

Delamination Onset in (O2/62^-02)s AS4/3501-6 
Graphite Epoxy Laminates 



6, degrees 

Pmax , N (lbs) 

1 5 

20 

25 

30 

Static Test 

1 .35 

1.57 

2.67 

5.25 


14.2 

9.7 

11.0 

20.7 

15,570 (3500) 

12.0 

10.3 

13.7 

23.3 

13,340 (3000) 

9.0 

14.0 

15.7 

32.0 


TABLE 3 - AS4/3501-6 Graphite Epoxy Lamina Properties 

En « 135 GPa (19.5 x 10^ Psi) 

E 2 2 = 11 GPa (1.6 x 10 6 Psi) 

G 12 = 5.8 GPa (0.847 x 10 6 Psi) 

vi 2 = 0.301 

ai = -0.41 xIO* 6 /°C (-0.23 x10- 6 /°F) 
a 2 = 26.8 xlO-6 /°C (14.9 xlO’6 /°F) 

Pi =0.00 /weight % 

P 2 = 5.56 xIO' 3 /weight % 
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TABLE 4 - Partial Local Delamination Geometries 


Configuration 


Fig. 26a 

_a_ 

sin0 

a tan0 

tan0 

Fig. 26b 

a cos0 

a cos0 sin© 

cos0 sin0 

Fig. 26c 

a[co S e--SML] 


sin0 (cos0 - -Sinfi_| 


' tan30 > 

■MHI 

* tan30/ 


TABLE 5 - Percentage Difference between 
G M and G M+T ‘ H for AH = 0.3% 



9, degrees 

P max , N (lbs) 

1 5 

20 

25 

30 


5.8% 

7.5% 

7.7% 

7.0% 


5.0% 

6.6% 

6.8% 

6.1% 

17,790 (4000) 

4.3% 

5.8% 

6.2% 

5.8% 

22,240 (5000) 

3.7% 

4.9% 

5.1% 

4.9% 


TABLE 6 - Range of Critical Normalized Delamination Lengths at the 
Free Edge of (O2/02/-02)s Graphite Epoxy Laminates 


Measured After Local 

Delamination Onset in Fatigue 

0, Degrees 

a/L 

1 5 

.034-. 106 

20 

.036-. 090 

25 

.029-. 063 

30 

.019-. 062 
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Fig. la. RADIOGRAPH OF (02/202/-202) s LAMINATE 
AFTER MATRIX CRACKING ONSET IN FATIGUE 
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Fig. 11. LOCAL DELAMINATION MODELS 
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FIG. 12 Normalized G for Uniform Local Delamination 
9/ -9 Interface of (0/0/ -e)„ Laminates 



> Degrees 



FIG. 13 Normalized G for Edge Delamination in 
0/ -0 Interface of (0/0/ -0) Laminates 
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Fig. 15. PARTIAL LOCAL DELAMINATION GEOMETRIES 
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FIG.16 Normalized G for Delamination in 15/-15 Interface 
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FIG. 17 Normalized G for Delamination in 20/-20 Interface 

of (0/20/-20)^ Laminate 
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FIG. 18 Normalized G for Delamination in 25/-25 Interface 

of (0/25/-25)^ Laminate 








FIG. 19 Normalized G for Delamination in 30/-30 Interface 

of (0/30/-30)^ Laminate 







FIG. 20 Normalized G for Local Delamination in 0 / -e Interface 

of (0/9/ -e) Laminates 
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FIG. 21 q M+T+H FOR LOCAL DELAMINATION 
FROM -15 DEGREE MATRIX CRACK 










FIG. 23 g m+t+h FOR LOCAL DELAMINATION 
FROM -25 DEGREE MATRIX CRACK 
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N, Cycles to Delamination Onset 
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